The WiggleZ Dark Energy Survey: 
Cosmological neutrino mass constraint from blue high-redshift galaxies 
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The absolute neutrino mass scale is currently unknown, but can be constrained by cosmology. 
The WiggleZ high redshift, star-forming, and blue galaxy sample offers a complementary dataset 
to previous surveys for performing these measurements, with potentially different systematics from 
non-linear structure formation, redshift-space distortions, and galaxy bias. We obtain a limit of 

m v < 0.60 eV (95% confidence) for WiggleZ+Wilkinson Microwave Anisotropy Probe. Combining 
with priors on the Hubble parameter and the baryon acoustic oscillation scale gives m„ < 0.29 eV, 
which is the strongest neutrino mass constraint derived from spectroscopic galaxy redshift surveys. 



I. INTRODUCTION 

Neutrinos are the lightest massive known particles, 
yet they are treated as exactly massless by the Stan- 
dard Model of particle physics. We know they have non- 
zero masses because neutrino oscillation experiments us- 
ing solar, atmospheric, and reactor neutrinos have mea- 
sured mass differences between the three species to be 
Am| 2 = 2.43 x HT 3 eV 2 and Amf^ = 7.59 x l(T 5 eV 2 
[l], The Heidelberg-Moscow experiment has limited 
the mass of the electron neutrino to be less than 0.35 eV 
using /?— spectroscopy [3j, but no current experiment has 
sufficient sensitivity to measure the absolute neutrino 
mass. 

Massive neutrinos affect the way large-scale cosmolog- 
ical structures form by suppressing the gravitational col- 
lapse of halos on scales smaller than the free-streaming 
length at the time the neutrinos become non-relativistic. 
This leads to a suppression of the small scales in the 
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galaxy power spectrum that we observe today, and con- 
sequently we can infer an upper limit on the sum of neu- 
trino masses from the matter distribution of the Universe 
4]. Combining with the lower limit of ^2 m v > 0.05eV 
provided by the mass differences from oscillation experi- 
ments, allows us to narrow the range of possible neutrino 
masses. 

The cosmic microwave background (CMB) provides an 
upper limit of m v < 1-3 eV [5, all limits are 95% con- 
fidence]. Combining with large-scale structure measure- 
ments such as the galaxy power spectrum B-lsll. ga laxy 
luminosity function [Ioj |. cluster mass function [llT |l2fl. 
or the scale of baryon acoustic oscillations [BAO, [5|, Q 
tightens the constraints to J2 m v ~ 0.3 eV by break- 
ing degeneracies with other parameters. Consequently 
neutrino mass constraints are important goals of current 
and future galaxy surveys such as e.g. Baryon Oscillation 
Spectroscopic Survey [13j, Dark Energy Survey [14J and 
Euclid [l5[. In this letter we use the galaxy power spec- 
trum from the WiggleZ Dark Energy Survey to constrain 
the sum of neutrino masses. 

The WiggleZ galaxy survey has several complementary 
aspects and potential advantages over previous surveys: 
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FIG. 1. The ratio between the simulated WiggleZ halo power 
spectrum (from GiggleZ) and the corresponding linear power 
spectrum (normalised to the first bin) at z = 0.2 (blue dashed) 
and z = 0.6 (black solid) and for luminous red galaxies (dotted 
red). The vertical lines indicate our fitting range of k — 0.02 — 
0.3 ftMpc" 1 . Non-linear corrections are clearly less significant 
for the high-redshift, low- bias WiggleZ halos than at lower 
redshifts. 



redshift range z < 1 [l7|- We split the data into four 
redshift bins of width Az = 0.2 with effective redshifts 
of z Q g = [0.22,0.41,0.6,0.78]. The power spectra, P ohs , 
and covariance matrices, C, are measured in Ak = 
0.01 /iMpc -1 bins using th e op timal-weighting scheme 
proposed by Feldman et al. [18| for a fiducial cosmolog- 
ical model [Tgj ] . The Gigaparsec WiggleZ Survey (Gig- 
gleZ) simulations [2(| were designed to probe the low- 
mass haloes traced by WiggleZ galaxies, whilst providing 
an equivalent survey volume allowing the measurement of 
power spectrum modes with k = 0.01 — 0.5 /iMpc" 1 and 
provide a powerful means for testing and calibrating our 
modeling algorithms. The GiggleZ simulations show that 
over the range of scales and halo masses relevant for this 
analysis, the galaxy bias is scale-independent to within 
1% [20j whereas the neutrino scale dependent effect is of 
the order of 5% for a 0.3 eV neutrino mass [2l| . 



1) The neutrino suppression of the galaxy power spec- 
trum is degenerate with effects from non-linear structure 
formation. Non-linearities increase with time so for the 
distant galaxies probed by WiggleZ, the contamination 
from non-linearities is smaller than for previous surveys. 
This is illustrated in Fig. [T] where we show the ratio be- 
tween a simulated WiggleZ power spectrum and the lin- 
ear power spectrum for z = 0.2 (dashed blue) and z = 0.6 
(solid black) . For comparison we also show the ratio for 
simulated highly-biased massive haloes at z = 0.2 (dot- 
ted red). 

2) The relationship (bias) between the observed galaxy 
distribution and the dark matter distribution, which is 
influenced by massive neutrinos, depends on the observed 
galaxy type. Previous studies [e.g. 0, HH measured red 
galaxies, which tend to cluster in the centers of dark mat- 
ter halos, whereas the star-forming blue WiggleZ galaxies 
avoid the densest regions. This leads to a lower overall 
bias, which makes WiggleZ less susceptible to any possi- 
ble systematics that could arise from a scale-dependence 
of the bias. 

Galaxy redshifts are not entirely determined by the 
Hubble flow, but also by their peculiar motions (redshift- 
space distortions) providing a challenge when comparing 
observations (redshift-space) with theory (real-space). 
Through exhaustive tests using numerical dark matter 
simulations of the WiggleZ survey, we demonstrate the 
breakdown of common models at small scales, and cali- 
brate a new non-linear fitting formula. 



II. THE WIGGLEZ DARK ENERGY SURVEY 

The WiggleZ Dark Energy Survey was designed to de- 
tect the BAO scale at higher redshifts than was possi- 
ble with previous datasets. The 238,000 galaxies are se- 
lected from optical galaxy surveys and ultraviolet imag- 
ing by the Galaxy Evolution Explorer to map seven re- 
gions of the sky with a total volume of 1 Gpc 3 in the 



III. METHOD 

Large-scale structure alone cannot determine all cos- 
mological parameters, so we include data from the CMB 
as measured by the Wilkinson Microwave Anisotropy 
Probe (WMAP). To compute the parameter likelihoods, 
we use importance sampling [22l . [23l j of the WMAP7 
Markov Chain Monte Carlo (MCMC) chains available 
online [Hj) from fitting to WMAP alone as well as 
to the chains combining WMAP with the BAO scale 
from SDSS Luminous Red Galaxies [l6| and a H = 
74.2 ± 3.6 km s -1 Mpc - prior on the Hubble parame- 
ter [25[ (BAO+7?o)- In order to disentangle the neu- 
trino mass and non-linear effects, we test six different 
approaches for generating a model redshift-space galaxy 
power spectrum. This section explains the aspects com- 
mon to all models. 

Matter power spectra: First we calculate the mat- 
ter power spectrum, P mi for each redshift bin for the set 
of cosmological parameters 9 = [J7 C (cold dark matter 
density), fib (baryon density), fl\ (dark energy density), 
il„ (neutrino density), h (Hubble parameter), n s (spec- 
tral index), A 2 R (amplitude of primordial density fluctu- 
ations)]. We assume a standard flat ACDM cosmology 
with no time variation of w in agreement with observa- 
tional data [H . The effective number of neutrinos is fixed 
to iV e ff = 3.04 assuming no sterile neutrinos or other rel- 
ativistic degrees of freedom (the 0.04 accounts for the 
non-thermal nature of the neutrino spectrum, which gets 
skewed during decoupling because higher-energy neutri- 
nos decouple later than lower-energy neutrinos) . The P m 
is converted to a galaxy power spectrum, -P ga i, using one 
of the six approaches described in Sec. IIVI 

Scaling and convolution: Before P ga i can be 
compared to the observed power spectrum, P b s , the 
survey geometry and the fiducial cosmological model 
used when measuring the power spectra must be 
accounted for by Alcock-Paczynski scaling, aj? cl = 
(D\H^ 1 (z))/(D\ &d H^(z)), and convolution with the 
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survey window function, Wij 23, 26, 27]: 

W i:j (k)P gal (kj/a scl ) 



satisfies: 



-^con(^i) ^ 



(1) 



where Z?^ is the angular diameter distance and H(z) is 
the Hubble parameter. Details of the window function 
can be found in (28|. For all models, we marginalise an- 
alytically [22l ] or numerically over a linear galaxy bias 
factor. 

Likelihood: We assume the power spectra to be dis- 
tributed as a multivariate Gaussian so the likelihood can 
be determined as: 



21n(L(0))= X 2 = ^A^ 1 A,, 



(2) 



where 9 is the set of cosmological parameters (including 
galaxy bias), A, = [P obs (fci, 6)-P con {ki, 9)] withP con (fci) 
being the convolved power spectrum in the i'th bin, and 
Cij the covariance matrix. 

The power spectrum measurements in the seven sur- 
vey regions are treated as independent observations, and 
their likelihoods combined by multiplication. We require 
the bias to be the same for all regions at a given red- 
shift. However, we allow the bias to vary between red- 
shifts since the survey magnitude and colour cuts cause 
the galaxy luminosities to evolve with redshift. 

Importance sampling: We use importance sampling 
to re- weight the WMAP MCMC likelihood chains [Hj]. 
For a chain of parameter values 9 drawn from a likeli- 
hood, L, it is possible to re- weight the likelihoods with an 
independent sample from the same underlying parameter 
distribution. The WMAP and WiggleZ power spectra are 
independent measurements, so their likelihoods can be 
combined by multiplication. Thus the weight, i^wmapj 
of each element in the MCMC chain, i, can be re- weighted 

by ^WMAP+WigglcZ = ^WiggleZ (#')<4vMAP [HE!- 

Using the CosmoMC software [29] for a subsample of 
the data, we have checked that the preferred regions of 
parameter space for WMAP and WiggleZ overlap, and 
consequently importance sampling is a valid method. A 
CosmoMC module for the WiggleZ power spectra is un- 
der development (30|. We have also fitted the WiggleZ 
power spectra alone over the range k = 0.02—0.2 h Mpc -1 
varying only Q m and where = f2i,/fi TO , keeping all 
other parameters fixed at the WMAP 7 best fit values 
i (W iggleZ alone can not constrain all the cosmological 
parameters). The resulting parameter values and un- 
certainties are consistent with the measurements of these 
parameters using WMAP data alone, which validates the 
assumption that the two data sets probe the same cos- 
mological parameter space. 

Neutrino mass constraint: The neutrino mass limit 
is calculated from the histogram of the WMAP MCMC 
chain likelihoods re- weighted by WiggleZ. The X% con- 
fidence upper limit on ^ m v is the value of M l j m that 
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Using the WMAP7 chains alone gives a 95% confidence 
limit on the neutrino mass of ^ m„ < 1.3 eV [3l| and 
< 0.55 eV when combining with BAO+Ho 



IV. MODELING APPROACHES 

Massive neutrinos suppress the power spectrum on 
all scales smaller than their free-streaming length at 
the time the neutrinos become non-relativistic. For 
^2 fn v = 0.3 eV the most significant suppression happens 
for k = 0.3 — 1.3/i Mpc -1 , but the fc-dependence of the 
suppression is most pronounced for k = 0.1 — 0.3 h Mpc -1 
and consequently easier to disentangle from other cos- 
mological parameters [1, [2lJ . At low redshift structure 
formation is no longer linear for k > 0.1 /iMpc -1 [il. I32I- 
|34| . The standard way of determining the matter power 
spectrum of non-linear structure formation is the phe- 
nomenological Halofit calculation [35j distributed with 
CAMB [36j. Halofit has been derived for massless neu- 
trinos, but using hydro-dynamical SPH simulations Bird 
et al. [2l| demonstrated that including realistic neutrinos 
in the simulations only produced a change in the power 
spectrum amplitude of less than 1% for k < 0.3 h Mpc -1 . 
Simulations show that redshift-space distortions become 
/c-dependent at low redshift and consequently are degen- 
erate with neutrino mass (37l . l38j . 

With the aim of constraining neutrino mass, Swan- 
son et al. [23| investigated 12 different models for non- 
linear structure formation, galaxy bias, and redshift- 
space distortions. They concluded that models with only 
one free parameter are unable to provide a good fit for 
U > 0.1 - 0.2 h Mpc" 1 for the SDSS red and blue 
galaxies. Blake et al. [28[ fitted 18 different models to 
the 2D redshift-space WiggleZ power spectrum for a fidu- 
cial cosmology, and concluded that the best fitting mod- 
els are those of Jennings et al. [37j and Saito et al. [39[ 
but the latter is very computationally expensive. With 
these conclusions in mind, we have tested six different ap- 
proaches for our analysis. The models are fairly similar 
at low values of k, where the large-scale clustering can be 
treated as linear. There the theory is quite robust, and 
we expect little difference between the models. However 
the difference between the models starts to increase for 
k > 0.2 /iMpc -1 , which significantly affects the outcome 
of the fitting. Throughout the analysis we have fixed the 
lower limit to be /c m j n = 0.02 h Mpc -1 , which corresponds 
to the largest modes observed in each of the WiggleZ re- 
gions (the final results are not very sensitive to the exact 
value). We then present all results as a function of fc max - 
The models are shown in Fig. [5] for a fixed cosmology 
and described below. Testing the models on simulations 
showed that the models A)-E) are insufficient and the 
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FIG. 2. A weighted average of the WiggleZ power spectra 
in the survey regions and redshifts, and the six models for 
the best fit cosmology of model F). The models are: A) blue 
dotted, B) green dashed, C) magenta dot-dashed, D) cyan 
triple dot-dashed, E) red long dashed, F) thick orange solid. 
The vertical lines are k = 0.02 h Mpc -1 and k = 0.3 h Mpc -1 . 
The divergence between the models at large k is clear and 
demonstrates the importance of careful modeling. 



complexity of model F) is necessary. 

A) Linear: We use CAMB to calculate the linear 
matter power spectrum, to which we add a linear bias 
model with redshift-space distortions in the Kaiser limit 
[40l |. The model is valid within a few percent for k < 

0.15/1 Mpc" 1 m. 

B) Non-linear: We use Halofit to calculate the 
non-linear matter power spectrum, Phf- The bias and 
redshift-space distortions are treated as for model A). 

C) Non-linear with fitting formula for redshift- 
space distortions and pairwise velocities: Combin- 
ing the ansatz of Scoccimarro [34j with fitting formulae 



derived from simulations, the model of Jennings et al. 37 1 
describes the time evolution of the power spectra of dark 
matter density fluctuations and their velocity divergence 
field. The details of our implementations of this model 
are given in Parkinson et al. (30j . 

D) Non-linear with fitting formula for redshift- 
space distortions and zero pairwise velocity 
damping: The fitting formulae used in model C) were 
derived for dark matter particles and not for halos. Set- 
ting all galaxy velocity dispersions to zero in model C) 
provides a better fit to the GiggleZ halo catalogues, and 
we have treated this special case as a separate model. 

E) Non-linear with pairwise galaxy velocity 
damping: Non-linear structure formation leads to in- 
creased peculiar galaxy velocities at low redshift, which 
damps the observed power spectrum. The effect can be 
described by the empirical model fill ]: 



P ga l(fc) = b 2 r P M (k) 



(i + £^ 2 ) 2 



i 



r dfj, 



(4) 



(kfa v fi) 2 

where / is the cosmic growth rate, fj, = k ■ z is the cosine 



of the angle between the wave vector, k, and the direction 
of the line of sight, z, and the one-dimensional velocity 
dispersion, a v (in units of Mpc), is given by |34j |: 



2 1 



3 (2tt) 



dk'P hn (k') 



(5) 



Setting a v = h^ 1 Mpc we recover model B). 

F) A^-body simulation calibrated approach: All 

the non-linear effects are present in an A^-body simula- 
tion for a fiducial cosmology and can be implemented 
following the approach of Reid et al. Iff] . For each trial 
cosmology: 



ptrial/iA i 2 ptrial 

-'gal \ K ) — u ^hf.nw 



(k)- 



P 



trial 



dampcd(^) ^GiggleZ 



«0 P$LJk) 



(6) 



where 



•TdampedW 



ap(fc)) 

(7) 

and /damp(fc) = exp(— (kcr v ) 2 ) with a v given by Eqn. [5] 
^GiggleZ (^) ^ s f° un( i from a 5 th order polynomial fit to the 
power spectrum of a set of halos in the GiggleZ simula- 
tions chosen to match the clustering amplitude of Wig- 
gleZ galaxies. P nw and Phf ,nw are the power spectra with- 
out the acoustic peaks, for the linear and Halofit power 
spectra respectively. They are calculated from a spline 
fit to the CAMB power spectra following the approach 
of Jennings et al. 3?| and Swanson et al. [23] . The fac- 
tor of b 2 in Eqn. [5] is related to galaxy bias. The second 
factor represents the smooth power spectrum of the trial 
cosmology. The third factor defines the acoustic peaks 
and their broadening caused by the bulk-flow motion of 
galaxies from their initial positions in the density field, 
and the fourth factor describes all additional non-linear 
effects in the A^-body simulation. 

Performance of the approaches: We tested the 
different approaches by fitting to the z = 0.6 power 
spectrum of a GiggleZ halo catalogue matching the clus- 
tering amplitude of WiggleZ galaxies to two sets of 2D 
parameter grids: Q m — ft, and £l m — n s , with the re- 
maining parameters fixed at the GiggleZ fiducial cos- 
mology values. We chose these grids because the pa- 
rameters are susceptible to degeneracies with neutrino 
mass. In both cases we obtain very similar conclusions, 
so here we only present the results of fi TO — ft- For 
fcmax < 0.2/iMpc _1 most of the models produce a good 
fit, whereas for /c max = 0.3 ft. Mpc -1 models B), C) and 
E) break down and give reduced \ 2 values above 1.5. 
The upper panel of Fig. [3] shows the x 2 f° r t ne fiducial 
GiggleZ cosmological parameters, which is a measure of 
how well the models recover the input parameters. The 
lower panel of Fig. [3] shows the difference between \ 2 °f 
the GiggleZ values and the best fit, indicating how far 
the best fit is from the input values. We assume that the 
A^-body simulation, which provides a complete census of 
the relevant non-linear effects, yields the most accurate 
clustering model. In this sense the good performance of 
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FIG. 3. Upper: Reduced \ 2 of models A)-F) fitted to the N- 
body simulation halo catalogue for the GiggleZ fiducial cos- 
mology values. In absence of systematic errors the models 
should recover the input cosmology with x 2 /dof = 1. Lower: 
Difference in reduced \ 2 values when using the GiggleZ fidu- 
cial cosmological parameters and the best fit values. The 
models are: A) blue dotted, B) green dashed, C) magenta 
dot-dashed, D) cyan triple dot-dashed, E) red long dashed, 
F) thick orange solid. 



FIG. 4. Upper: Reduced \ 2 as a function of fe max for each 
of the six approaches. Lower: Upper limits on Yl m v as a 
function of fc max . The models are: A) blue dotted, B) green 
dashed, C) magenta dot-dashed, D) cyan triple dot-dashed, 
E) red long dashed, F) thick orange solid.. The dashed grey 
line is the lower limit from oscillation experiments, and the 
black lines are upper limits from WMAP+BAO+Ho (dotted) 
and WiggleZ+WMAP+BAO+Ho (solid). 



model F) (Fig. [3]) is a consistency check, and the vari- 
ations of results produced by the other models are due 
to the breakdown in their performances compared to the 
simulation. We are cautious about fitting too small scales 
where modeling of the non-linearities and massive neu- 
trinos become less robust, and the data is dominated by 
shot noise. The neutrino implementation in CAMB is 
only accurate for fc max < 0.3 ft Mpc" 1 [21] which we take 
as an upper limit for our fits. 



V. RESULTS AND DISCUSSION 

When fitting the observed WiggleZ power spectra we 
obtain the results presented in Fig. |4j The upper panel 
shows the \ 2 as a function of fc max for the best fit- 
ting parameter values for each of the six approaches, 
and the lower panel shows the corresponding neutrino 
mass constraints. Although all models produce sim- 
ilar x 2 values, our comparison with the full A^-body 
simulation catalogue (Fig. [3]) revealed that systematic 
errors arise when models A) to E) are fit across the 
range of scales fc max < 0.3 ft. Mpc -1 . Using the fully- 
calibrated model F), we obtain < 0.60 eV for 
WMAP+WiggleZ with /c max = 0.3 h Mpc" 1 . Combining 
with BAO+-£/o reduces the uncertainty in il m and Ho, 
leading to stronger neutrino mass constraints. Without 
WiggleZ, the WMAP+iJo+BAO dataset gives £>„ < 
0.55 eV whereas combining with WiggleZ adds informa- 
tion about the power spectrum tilt (n s ). The resulting 



neutrino mass constraint is X) TO f < 0.29 eV for model 
F) and fc max = 0.3 ft Mpc" 1 . 

The relative probability distributions of m„ for 
model F) with fc max = 0.3 ft Mpc -1 are shown in Fig. [5] It 
is clear how adding WiggleZ data to the fit narrows the 
distributions (dotted to solid) both with (orange) and 
without (black) the inclusion of BAO+ifo- This is the 
strongest neutrino mass limit so far derived from spectro- 
scopic redshift galaxy surveys. The advantages of Wig- 
gleZ are a higher redshift for which the structure forma- 
tion is linear to smaller scales, and a simple galaxy bias 
for the strongly star-forming blue emission line galaxies. 

Our result is comparable to that obtained using pho- 
tometric redshift galaxy surveys [J2 m i' < 0.28 eV, 0,0], 
but the systematics in the two data set are completely 
different. For example, imaging surveys are potentially 
susceptible to systematic errors from the imprint of stars 
on the selection function [42j and the shape of the redshift 
distribution. WiggleZ contains negligible star contami- 
nation and much higher redshift resolution compared to 
photometric surveys. The high redshift and blue galax- 
ies of WiggleZ allow us to fit the power spectrum to 
smaller scales than previous surveys (both spectroscopic 
and photometric), where the effect of the neutrinos is 
larger, and get a similar neutrino mass constraint from 
a smaller, but well understood galaxy sample. Also the 
result from galaxy clusters {J2 m v < 0.33 eV, [ill, ^ s 
similar, but with different systematics. 

Since the data sets are all independent, they can poten- 
tially be combined in the future to provide even stronger 
constraints. This is particularly interesting in light of re- 
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2m v [eV] 

FIG. 5. The relative probability distribution of V} m„ 
from fitting model F) with fc max = 0.3/iMpc -1 for 
WMAP (dotted orange), WiggleZ+WMAP (solid 
orange), WMAP+BAO+#o (dotted black) and 
WiggleZ+WMAP+BAO+tfo (solid black). The dashed 
grey line is the lower limit from oscillation experiments, and 
the vertical lines are 95% confidence upper limits. 



cent results [a, |43j that hint at the existence of additional 
neutrino species (N c g > 3.04). Allowing for additional 
neutrino species degrades the constraining power of large 
scale clustering alone, and the combination of N e g and 
^2 Tn-u is therefore poorly constrained with current data. 

In the future, galaxy surveys such as the Baryon 
Acoustic Oscillation Survey, Dark Energy Survey and 
Euclid will be far more sensitive, giving cosmological neu- 



trino mass constraints of order J2 m v < 0-05 — 0.1 eV 
0, HU- This will be small enough to distinguish be- 
tween the ordering of the neutrino masses (normal hi- 
erarchy where mi < « or inverted where 
77i 3 << mi < "12)- However, as demonstrated in this 
paper, the small details of the modelling of non-linear ef- 
fects become very important, so robust modelling either 
theoretically or calibrated to simulations with massive 
neutrinos will be necessary. 
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